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A numerical  procedure  is  used  to  examine  the  influence  of 
transverse  shear  forces  in  the  yield  criterion  and  rotatory 
inertia  on  the  dynamic  plastic  response  of  beams.  Various 
results  are  presented  for  a long  beam  impacted  by  a mass  and 
a simply  supported  beam  loaded  impulsively,  both  of  which  are 
made  from  a rigid  perfectly  plastic  material  with  yielding 
controlled  by  the  Ilyushin-Shapiro  yield  criterion. 

Transverse  shear  effects  lead  to  a dramatic  reduction  in 
the  slopes  of  the  deformed  profiles  for  both  beam  problems. 
Moreover,  the  slope  of  the  deformed  profile  underneath  the 
striker  in  the  impact  problem  is  quite  sensitive  to  the  actual 
shape  of  a yield  curve,  while  the  maximum  transverse  displace- 
ment is  less  sensitive.  The  retention  of  rotatory  inertia  in 
the  basic  equations  leads  to  further  reductions  up  to  17  and 
10  per  cent  in  the  slopes  and  maximum  transverse  displacements, 
respectively. 


ii 


Notation 

(a)  Notation  used  throughout  the  article 

e,  f defined  by  equations  (6)  and  ( 17 f ) , respectively, 
m mass  per  unit  length  of  a beam 
t time 

w transverse  displacement  defined  in  Figure  1. 

x axial  coordinate  defined  in  Figures  1,  4(b),  and  6(b). 

z location  of  interface  in  equations  (2) 

H beam  thickness 

2 

Ir  mk  , where  k is  radius  of  gyration 

M,Q  bending  moment  and  transverse  shear  force  (per  unit 
length)  defined  in  Figure  1 

Mq  fully  plastic  bending  moment  capacity  of  cross-section 

Qo  fully  plastic  transverse  shear  capacity  of  cross-section 

M,Q  M/Mq,  Q/Qo 

Y shear  angle 

k curvature 

rotation  of  mid-plane  due  to  bending 
ft  defined  by  equation  (9b) 

(’)  3/3t , 3/3T. 


(b)  Additional  Notation  used  in  Section  4 

z location  of  interface  defined  in  Figure  4(b). 

G mass  of  striker 

I2  mIr/G2 

defined  by  equation  (17g) 
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2 

R energy  divided  by  GVq  /2 

T dimensionless  time  defined  by  equation  (]7i) 

dimensionless  duration  of  first  phase  of  motion 
initial  velocity  of  mass  G 

I 

V , v,  V2  transverse  velocities  of  beam  defined  in 

Figure  4(b). 

v',  v,  V2  v'/vo_  v/vQ,  v2/v0. 

W 12mH  W/G2v  2. 

o o 

• 

V 2 defined  by  equation  ( 1 7 j ) 

a x/z 

0 angular  rotation  at  centre  of  one  half  of  beam 

e 6m  6/gv  2 

o o 

\,  \Q  defined  by  equations  (17b,  c) 


(c) 


K 

2L 
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Additional  Notation  used  in  Section  5 

locations  of  interfaces  in  Figure  6(b). 

lr/mL2 

defined  by  equation  (54d) 
span  of  beam 

2 

bending  energy  at  time  t divided  by  mLVQ 

2 

kinetic  energy  at  time  t divided  by  mLVo 

energy  absorbed  in  shear,  g deformations  at  time  t 
2 

divided  by  n\LVQ 

dimensionless  time  defined  by  equation  ( 54 j ) 
dimensionless  duration  of  first  phase 
initial  impulsive  velocity 
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V,  defined  in  Figure  6(b) 

v,  vx  v/vQ,  V;L/vo 

W M w/mL2V  2 

o o 

WT  maximum  permanent  dimensionless  transverse  displacement 

a x/z^ 

®1 ' % ^/h , z 2/L 

6 angular  rotation  at  supports 

"6  M 0/mLV  2 

o o 

v Q L/2M 

o o 
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1.  Introduction 

It  was  remarked  in  References  [l]  and  [2]  that  the 
importance  of  transverse  shear  effects  has  recently  become  evident  in 
the  dynamic  plastic  response  of  various  structures.  However, 
the  available  theoretical  solutions  on  the  dynamic  plastic  behavior 
of  beams  which  retain  the  influence  of  transverse  shear  employ  a 
square  yield  condition  relating  the  bending  moment  (M)  and  trans- 
verse shear  force  (Q)  required  for  plastic  flow  [1,J,4J.  jt  is  the 
objective  of  this  article  to  examine  the  influence  of  transverse 
shear  and  rotatory  inertia  on  the  dynamic  plastic  response  of  1 earns, 
the  behavior  of  which  is  controlled  by  more  realistic  yield  cri- 
teria relating  the  transverse  shear  force  and  bending  moment. 

A simple  numerical  procedure  which  employs  piecewise 
linear  approximations  with  any  desired  accuracy  for  any  yield 
criterion  is  developed  herein,  and  used  to  reconsider  the 
influence  of  transverse  shear  and  rotatory  inertia  effects  on 
the  two  dynamic  beam  problems  investigated  in  Reference  [ 1 ] . 

2.  Basic  Equations 

The  equilibrium  equations  for  a Timoshenko  beam  may  be 
written  [1,3] 

3Q/3x  = mw  - q and  Q + 3M/3x  = Ir4i  (la,b) 

where  the  total  slope  of  the  centre  line  of  a beam  is 
3w/3x  = ip  + Yf  ^ is  the  rotation  of  the  line  elements  along 
the  centre  line  due  to  bending,  y is  the  shear  angle  of  points 

1 

I along  the  centre  line,  and  the  remaining  symbols  are  defined 

in  Figure  1 and  the  accompanying  Notation. 
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It  may  be  shown  [ 1 , 3 ] that  continuity  of  the  transverse 
displacement  (w)  and  angular  deformation  (ij>)  associated  with  a 
moving  interface  in  a rigid  perfectly  plastic  Timoshenko  beam 
requires  [1,3] 

[Q]  = mz2 [y]  (2a) 

and  [M]  = Irz2 [k]  , (2b) 

where  k = dip/dx,  and  [X]  means  the  difference  of  X on  either 
side  of  an  interface  which  travels  with  a velocity  z.  Various 
other  kinematic  and  equilibrium  requirements  at  travelling 
interfaces  are  given  by  equations  (3)  to  (5)  in  Reference  [1] . 

There  has  been  considerable  discussion  in  the  literature 
and  no  clear  resolution  as  to  whether  an  interaction  curve 
relating  the  generalised  stresses  M and  Q is  a proper 
yield  criterion  or  not.  This  question  has  been  reconsidered 
in  Reference  [5]  by  examining  all  previous  objections,  including 
Heyman's  [6],  and  it  was  concluded  that  a yield  curve  could 
be  used  within  the  framework  of  beam  theory.  It  was 
suggested  that  a suitable  compromise  between  the  simple  local 
(stress  resultant)  and  more  rigorous  non-local  (plane  stress, 
plane  strain)  theories  might  be  achieved  for  I-beams  when  using 
a local  theory  with  a maximum  transverse  shear  force  based  only 
on  the  web  area.  If  this  viewpoint  is  accepted,  then  it  is  not 
possible  to  describe  the  local  behavior  of  beams  at  supports, 
or  underneath  concentrated  loads,  which  is  similar  to  the 
prevailing  situation  with  classical  elastic  beam  theory. 

It  is  evident  from  Figure  5 in  Reference  [5]  that  a 


number  of  local  and  non-local  theories  predict  similar  curves 
in  the  M/Mq  - Q/Qq  plane  for  beams  with  rectangular  cross-sections. 
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Ilyushin  [7]  developed  a yield  surface  for  thin  shells 
with  solid  cross-sections  using  the  von  Mises  yield  criterion 
and  the  usual  assumptions  of  thin  shell  theory.  Shapiro  [8] 
extended  this  theoretical  work  in  order  to  cater  for  transverse 
shear  effects.  It  might  be  demonstrated  when  using  standard 
relations  between  inverse  hyperbolic  and  logarithmic  functions 
that  Hodge's  [9]  yield  curve  for  a beam  with  a rectangular  cross- 
section  is  identical  to  Ilyushin-Shapiro ' s theoretical  predictions 
[7,8,10]  for  the  plastic  yielding  of  a beam  which  is  subjected 
to  a bending  moment  and  a transverse  shear  force.  This  yield 
criterion  is  drawn  in  Figure  2 and  may  be  written 

Q/Qq  = <J>/sinh<|>  (3a) 

-2 

and  M/Mq  = coth<j>  - <$>sinh  (J>  (3b) 

where  4>  is  a parameter. 

The  theoretical  procedure  developed  herein  may  be  used 
for  any  yield  criterion,  but  the  Ilyushin-Shapiro  yield  criterion 
(equations ( 3) ) is  regarded  as  exact  in  this  article. 

3.  Details  of  the  Numerical  Procedure 

The  yield  criterion  in  the  proposed  numerical  procedure 
is  replaced  by  an  inscribing  or  circumscribing  polygon  having 
"n"  sides,  or  facets,  as  illustrated  in  Figure  2 for  the 
Ilyushin-Shapiro  yield  criterion  governed  by  equations  (3)  . 

These  piecewise  linear  approximations  coincide  with  the  exact 

yield  curve  when  n -*■  °°  . It  is  evident  from  the  corollaries  of 
the  limit  theorems  of  plasticity  that  the  static  collapse  loads 
associated  with  inscribing  and  circumscribing  yield  criteria 


- A 
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are,  respectively,  lower  and  upper  bounds  to  the  static  collapse 
load  according  to  the  exact  yield  condition.  The  choice  of  n 
is  therefore  related  to  the  acceptable  difference  between  the 
two  bounds . 

Now  it  is  apparent  from  Figure  2 that  two  types  of  plastic 
deformation  are  associated  with  a piecewise  linear  yield  condition: 
Type  A.  The  transverse  shear  force  (Q^)  and  bending 
moment  (M^)  are  known  when  the  generalised  stresses  are 
located  at  a node  "i" . However,  the  generalised  plastic 
strain  rate  vector  is  unknown  but  must  lie  within  the 
normals  associated  with  the  two  adjacent  facets  "i  - 1" 
and  "i". 

Type  B.  If  the  generalised  stresses  lie  on  a facet  "i", 
then  the  generalised  plastic  strain  rate  vector  must  be 
normal  to  this  facet.  The  generalised  stresses  can  lie 
anywhere  on  the  facet  "i"  between  the  two  adjacent  nodes 
"i"  and  "i  + 1”. 

A secant  linearisation  of  the  exact  yield  curve  is  now 
constructed  for  the  straight  line  which  cuts  the  exact  yield 
condition  as  shown  in  Figure  3.  The  equation  of  this  line  is 

0/Qo  + (bs/as> (M/Mc)  - bs  , (4) 

which  may  be  written 

Q/Qo  + <bs/ag) (M/Mo)  = i bg  (5) 

in  order  to  describe  a straight  line  which  cuts  the  yield 
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curve  in  any  of  the  four  quadrants,  where  ag  and  bg  are  the 

absolute  values  of  the  intersections  with  the  M/M  and  Q/Q 

o o 

axes,  respectively. 

A tangent  linearisation  of  the  exact  yield  curve  is 
constructed  using  straight  lines  which  are  tangential  to  the 
exact  yield  criterion  as  indicated  in  Figure  3. 

4.  Impact  of  a Mass  on  a Long  Beam 
4.1  Introduction 

Consider  the  infinitely  long  rigid  perfectly  plastic 
beam  shown  in  Figure  4a  which  is  hit  by  a mass  G travelling 
with  an  initial  velocity  V . The  influence  of  rotatory 
inertia  and  transverse  shear  forces  on  the  dynamic  plastic 
response  of  this  particular  beam  was  examined  from  an  analytical 
viewpoint  in  Reference  [1]  with  the  aid  of  a square  yield 
condition  relating  Q and  M.  A numerical  solution  is  sought 
in  this  section  using  the  Ilyushin-Shapiro  yield  curve,  which 
is  given  by  equations  (3)  and  illustrated  in  Figure  2.  The 
response  of  this  beam  consists  of  a number  of  phases  which  are 
described  in  the  following  sections. 

4.2  First  Phase  of  Motion 

Now,  guided  by  the  analytical  solution  presented  in 
Reference  [1] , it  is  conjectured  that  the  velocity  profile  shown 
in  Figure  4b  governs  the  behavior  of  this  problem.  A stationary 
hinge  develops  underneath  the  striker  at  x = 0,  the  region 
0 5.  x S.  z rotates  as  a rigid  body,  while  z S.  -x.  < z + i is  a 
plastic  zone,  and  x > z + i is  rigid  and  stationary. 

The  theoretical  results  in  Reference  [1] 


6 


indicate  that  the  bending  moment  is  MQ  and  Q is  negative 

and  quite  small  in  the  plastic  zone  z <.  x <.  z + l.  This 

suggests  that  the  plastic  zone  might  be  confined  to  the  vicinity 

of  the  M/Mq  axis  of  the  yield  curve  in  Figure  2 and  that  the 

piecewise  linear  yield  condition  shown  in  Figure  5 might  be 

adequate.  Equation  (5)  therefore  gives 

M/Mq  = 1 + eQ/Qo  (6) 

since  a =1  and  where  e = 1/b  . If  e = 0,  then  the  tangent 
s s 

line  M/Mq  = 1,  which  was  used  in  Reference  [1],  is  recovered, 
while  e > 0 corresponds  to  an  inscribed  facet.  These  two 
cases  would  bound  the  predictions  of  the  exact  yield  surface 
quite  accurately  when  Q/Qq  remains  small. 

Normality  of  the  generalised  strain  rate  vector  to  the 
yield  condition  represented  by  equation  (6)  requires 

Y/k  = -eK0/Qo  (7) 

• • • 

which  may  be  used  to  express  3w/8x  = ^ + y in  the  form 

9w/9x  = ip  - e(Mo/Qo)  9^/9x  . (8) 

It  may  be  shown  that  equation  (la)  with  q = 0 
and  equations  (lb) , (6)  and  (8)  give  the  partial  differential 

equation 


92i()/9x2  - R = 0 , 

(9a) 

where 

2 r,  2 , 2 2 _ _ 2.  -1 

Q = mQ0  (me  Mq  + IrQQ  ) 

(9b) 

Thus , 

equations  (9a)  and  (8)  predict 

i/j  = B sinhftx  + Ccoshftx 

(10) 

and 

w = fi-1(B  - eMonC/QQ)  coshnx  + 

+ f<_1(C  - eMQftB/Qo)  sinhQx  + D , 

(ID 

where  B(t),  C(t),  and  D(t)  remain  to  be  evaluated  from  the 
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boundary  conditions  and  continuity  requirements.  If  Q and 
M are  evaluated  from  equations  (la) , (lb) , (10)  and  (11)  and 

substituted  into  equation  (6) , then  it  is  found  that 

D = 0 . (12) 

Now  integration  of  equations  (10)  and  (11)  with  respect  to 
time  gives 

ip(x,t)  = E sinhftx  + Fcoshftx  (13a) 

and  w(x,t)  = n_1(E  - eMQf)F/Qo)  coshfix  + ft-1  (F  - eMQflE/Qo)  sinhfix 

(13b) 

when  satisfying  the  initial  conditions  ip(xr0)  = w(x,0)  = 0 , 
and  where  E(t)  and  F(t)  are  unknown.  Equations  (13)  and 
the  requirements  w(z,t)  = V2(t)  , w(z  + L,t)  =0  and 

i(p(z  + t,  t)  = 0 give 

w(x,t)  = V2  exp.  j-ft(x  - z )|  (14a) 

and  i£(x,t)  = - (mA^/I^n)  (1  - eMo£2/QQ)  exp.  { -ft  (x  - z)\  , (14b) 

when  using  equation  (9b)  and  where  z <.  x <.  z + £ and  £-*<». 
Thus,  equations  (la) , (lb) , (14a)  and  (14b)  may  be  solved 
to  give  the  dimensionless  transverse  shear  force  and  bending 
moment  in  the  form 

Q (a  ,T)  = -2  A K_1  V.exp.  {KA  (1  - a) }-2A  V-A  exp . {KA ( 1 - a) } 

o z o z (15) 

and  M(a,T)  = 1 - 12fK_^  V2  exp.  {KA  (1  - a)  }-12V2AfK  exp.  {KA  (1  - a)} 

(16) 

when  1 £.  a £ 00  and  where  the  two  "constants”  of  integration 
are  found  by  requiring  that  the  yield  condition  (6)  is  satisfied 
throughout  the  plastic  zone  (1  < a < »)  at  all  times,  and 
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a = x/z,  X = mz/G  , XQ  = 6mMQ/QoG  , Q = Q/Qq  , M = H/Kq  , 

f = eMQft/Qo  , K = ftG/m  , V2  = V2/Vq  , T = 12mMot/G2Vo  , 

— * 2 • • 

V.,  = V,G  /12mM  , and  X = GV  z/12M  . (17a-k) 

z 2 o o o 

The  velocity  field  in  the  rigid  region  0 £ x £ z 

is  linear  and  may  be  written 

w(x,t)  = (mV2/Irn)  {Ir&ym  + (1  - eMQn/Qo)  (z  - x)  } (18) 

when  w(z,t)  = V2  and  noting  that  9w(z,t)/3x  in  the  region 
• + • 

0 £ x £ z is  equal  to  i|;(z,t)  given  by  equation  (14b).  Thus, 
equations  (la) , (lb)  and  (18)  give 

Q (a,T)  = 2Xq^2  {(1  - f) (1  - a/2) aX2  + I2KXa}/I2K  + 

+ 2XqV2XX(1  - f ) a/I 2K  + Q ( 0 ,T)  (19) 

and  M(a,T)  = -12V2  {(1  - f)  (1  - a/3)X3a2/2  +I2Ka2X2/2  }/I2  K - 

- 6V2X  a2X2(l-f)/I^<  - 12V2(1  - f)aX/K  - 

- 6aXQ(0,T)/X  ■+  M(0,T)  , (20) 

t If  [9w(z,t)/9t]  = 0 , then  equation  (5a) 

of  Reference  [1]  gives  [Q(z,t)]  = 0. 

Now  [Q(z,t)]  = 0 implies  [M(z,t)l  = 0 for  a non-linear  yield 

curve  when  disregarding  the  possibility  that  M has  opposite  signs 

on  either  side  of  an  interface.  Thus,  equation  (5b)  of  [1] 

• • • 
gives  [i//(z,t)J  = 0.  Finally,  if/(x,t)  = 9w(x,t)/3x  throughout 

the  rigid  region  since  y(x,t)  = 0,  from  which  follows  the 


continuity  requirement  used  to  obtain  equation  (18) . 
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and  V = 1 - T/A  , (28) 

o 

where  2 (1  - f)  + |6I2K  - (1  - + 2I2  ( 6f -KAQ)  ^ + 

+ 12fI2/K  - 2I2Xq  = 0 , (29) 

and  V = w(0,t)/VQ 

Equation  (5)  can  be  written  Q/Qq  + (b1/a^)(M/MQ)  = -b^ 

for  the  plastic  hinge  immediately  underneath  the  striker. 

Thus,  normality  of  the  generalised  strain  rate  vector  to  this 
facet  demands 

k/Y  = (bj/a^  (Q0/Mq)  (30) 

which  becomes 

(V  - V2)/(V  - V)  = (6X1/Xq) (b1/a1)  (31) 

since  y = - (V  - V)/£  and  k = - (V  - V2)/z l for  a plastic 
hinge  of  length  "l". 

The  first  phase  of  motion  continues  at  the  node 


M ( 0 , T)  = 0,  Q ( 0 ,T)  = -1  until  the  dimensionless  time  T^  when 
the  generalised  strain  rate  vector  in  the  beam  underneath  the 
striker  becomes  normal  to  the  adjacent  facet  of  the  piecewise 
linear  yield  surface.  If  equations  (26)  to  (28)  are  substituted 
into  the  left  hand  side  of  equation  (31) , then  the  dimensionless 
time  T^  occurs  when  equation  (31)  is  satisfied,  or 

T1  = X 01(1  " f,Al2  + 2IS  + 2l2({  (1  ~ f)X0(a1/b1)/6  + 

+ (1  - f)  (1  + A.)\.  +(2A.  + l)i2K+2[2r1  . 

A 1 1 (32) 


A second  phase  of  motion  commences  at  the  dimensionless 
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time  with  the  generalized  stresses  at  the  plastic  hinge 
underneath  the  striker  lying  on  facet  1 in  Figure  5.  This 
phase  of  motion  with  X^o  is  eventually  completed  at  T2  when 
the  generalized  stresses  have  migrated  along  facet  1 and  reached 
node  2.  A third  phase  of  motion  then  commences  at  T2  with  the 
generalized  stresses  remaining  at  node  2 and  the  generalized 
strain  rate  vector  rotating  with  time  until  it  becomes  normal 
to  the  adjacent  facet  2 at  a dimensionless  time  • A fourth 
stage  of  motion  then  follows  with  the  generalized  stresses 
migrating  along  facet  2 in  Figure  5.  It  is  evident  that  the 
response  of  this  beam  consists  of  a sequence  of  phases  which 
are  characterized  as  either  type  A or  type  B which  are  defined 
in  section  3. 

The  governing  equations  are  sought  in  the  next  two 
sections  for  the  cases  when  the  generalized  stresses  are 
either  at  a typical  node  i (type  A)  or  on  a typical  facet  i 
(type  B)  of  the  piecewise  linear  yield  condition  in  Figure  2. 

4.3  Type  A Response 

As  remarked  earlier,  the  generalized  stresses  remain  at 
a node  during  a response  characterized  as  type  A.  Now,  Q(0,T)  = 

-d ^ and  M(0,T)  - -g ^ at  the  typical  node  " i"  in  Figure  2.  Thus, 
equations  (22)  and  (23)  may  be  integrated  to  give 

V2|l-f)  A2/I2K  +2  A + 2/K  | = dL  T/Ao  + (33a) 

and  2V2  |(l-f)  A3A2K  + 3A2  + 6f  A/K-*€f/K2  } = (1  + gA)T  + C2,  (33b) 
where  and  C2  remain  to  be  found.  If  the  present  phase  of 
motion  commenced  at  a dimensionless  time  T^  with  associated 
values  A^  and  v2(k)'  fc^en 

I 


(34a) 


C1  " v2(k)|(1*f>V  K+  2\  + 2/K)  ~ diVXo 

and  C2  = 2V2(k){(l-f)X^/I2K  + 3A2  + 6fAk/K  + 6f/K2}  - 

- (1  + gi)Tk.  (34b) 

Furthermore,  eliminating  V2  from  equations  (33)  gives  X at  a 
dimensionless  time  T in  the  form 

{2(l-f)(di  T/X0  + C1)/I2K^X3  + |6(d.T/Xo  + Cx)  - (1-f)  (1  + gi)T/l2K 
-(l-f)C2/I2K|x2  + 2|6f(diT/X0  + C-^/K  -(1  + gi ) T - C2|  X 
+ 12f  (d.T/X0+C1)/K2  - 2 1 ( 1 + g.  )T  + C2|/K  = 0.  (35) 

Thus,  equation  (35)  with  C-^  and  C2  given  by  equations  (34) 
predicts  X at  any  T which  then  allows  V2  to  be  calculated  using 
equations  (33a)  and  (34a) . V may  then  be  evaluated  from  equation 
(18)  with  x = 0,  or 

V = V2  I (1-f) A + I2k|/I2K  (36) 

* • 

while  V is  given  by  equation  (25) . Finally,  V2  and  X 

follow  by  differentiating  equations  (33) , or 

X = [<li/Xo  - ^2  j(l-f)A2/i2K  + 2 A + 2/k}]T2V2  {(l-fJX/^K  + l}]  , (37a) 

and  i/2  = [(1  + gi  )Ao{(l-f)AA2K  + l}  -3d.  {(l-f)X2/I2K  + 

+2X  +2f /K  |]  [ 2Xq|  1 + (1-f)  Vl2K|{(l-f)X3/]2K  + 

+3A2  + 6f X/K  + 6f/K2}  - 3Xq  {(1-f)  X2/32K  + 

+2X  + 2f/K  | {(1-f) A 2/I2  K + 2X  + 2/K  } ] . (37b) 

' \ 

This  phase  of  motion  commences  at  a dimensionless  time 
Tk  when  the  generalized  strain  rate  vector  underneath  the  striker 
is  at  node  i and  is  assumed  to  be  perpendicular  ^ to  facet  i-1 
of  the  piecewise  linear  yield  curve  in  Figure  2.  The  general- 
ized strain  rate  vector  remains  at  node  i and  rotates  during 
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the  ensuing  motion  until  the  phase  is  completed  at  a dimension- 
less time  T^+^,  when  the  generalized  strain  rate  vector  is 
perpendicular  to  facet  i or  returns  to  its  initial  position, 
which  is  perpendicular  to  facet  i-1.  Thus,  equation  (31)  gives 


‘vk  - v2(k)'lvk  - V <6VV 


(38a) 

J.  A X 

when  the  phase  commences  at  T^,  and  either 

I 

(V(k+l)  ' V2(k+l))/(V(k+l)  " V(k+1)>  = (6X(k+1)/X0)  (bi/ai), 
or 


(V(k+l)  * V2(k+l))/(V(k+l)  " V(k+1))  = (6x(k+i)/V  <bi-i/ai. 

when  this  phase  of  motion  is  completed  at  Tk+^.  V2(k4l)  an<* 

A (k+i ) are  the  initial  values  for  the  following  phase  of  motion 
which  commences  at  Tk+^. 


* 2.  Tk+i*  This  is  added  to  the  transverse  displacement 


The  transverse  velocities  within  the  two  zones  c <_  x < z 
and  z < x <_  °°  given  by  equations  (18)  and  (14a)  may  be  integrated 
numerically  with  respect  to  time  to  predict  the  contribution  to 
the  displacement  profile  accumulated  during  the  phase  of  motion 
Tk  < T < 

profile  at  T^  to  give  the  total  transverse  displacement. 

The  bending  moment  and  transverse  shear  force  distribu- 
tions in  both  zones  at  a dimensionless  time  T are  calculated 
from  equations  (15),  (16),  (19),  (20),  (33a),  (34a),  (35),  (37a) 
and  (37b) . 


If 


0,  T =0,  and  V_ , . = 0 when  motion  commences, 
o 2 (o) 


then  Tk+1  = Tjy  Ak+1  = Aj^,  and  = V2(i)  are  associated 

with  the  first  phase  of  motion  considered  in  section  4.2  provided 


i - 1,  di  = 1,  gt  - 0,  A = 0,  and  V i 


is  identified  with  V~. 


(38b) 
x) (38c) 
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It  is  straightforward  to  show  that  = C2  = 0 according  to 
equations  (34a)  and  (34b) , while  equations  (33a) , (35) , (36) , 
and  (38b)  reduce  to  equations  (26) , (29) , (27)  , and  (31)  , 
respectively.  Thus,  equation  (32)  for  is  recovered  since 

I 

V is  given  by  equation  (25)  . It  is  not  necessary  to  satisfy 
equation  (38a)  at  Tq  = 0. 

4 . 4 Type  B Response 

As  remarked  earlier,  the  generalized  stresses  lie  on  a 
facet  of  a piecewise  linear  yield  curve  during  a type  B response. 
Consider  the  facet  i in  Figure  2 for  which  equation (5)  may  be 
written 

Q(0,T)  + ( b ^a  i) M ( O , T ) = -bi  (39) 

for  the  plastic  flow  in  the  beam  directly  underneath  the  striker. 

The  associated  normality  requirement  according  to  equations  (30) , 

(31),  and  (38b)  is 

I 

(V  - V2)/(V  - V)  = (6A/Aq) (bi/ai) . (40) 

Now,  substituting  equations  (25)  and  (36)  into  equation 
(40)  gives 

V2  = I2K|(ai/bi)  (1-f)  Aq/6  + (1-f ) (1+A)  A +I2K(1+2A)  + 2I2}  (41) 

and  therefore  according  to  equation  (36) 

V=  {(1-f)  A + i2K}|(a./bi)  (1-f)  Ao/6  + (1-f ) (1+A)  A + I2K  (1+2A)  +2I2}  .(42) 

Equations  (22) , (23)  and  (39)  yield 

T = V2  [2  (1-f)  A3/I2K  +6A2  + 12f  A/K  + 12f/K2  + (a./b^  Aq  { (1-f ) A2/I2K  + 

+ 2A  + 2/K  } 1(1  + ai)  - Cx  (1  + a-^  , (43) 


where  Cx  = V2(k)  [ 2 (1-f ) Ak/I2K+  6A2  -t  12fAk/K  + 12f/K2  + 
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+ (a./b.  )A0  {(l-f)X£/i2K  + 2Ak  + 2/K}  ]-  (1  + a^T^.  (44) 

is  obtained  from  the  requirement  that  V2  = V2 ^ and  X = X^ 
at  T = T^  when  this  phase  of  motion  commences.  The  time 
derivatives  of  equations  (41)  and  (43)  give 

V2  =-  I2KX  (2  ( 1 — f ) X + 1-f  +2l2K}  { ( 1 — f ) X2  + (1-f  + 2l2K)X  + 

+ (l-f)X0(ai/bi)/6  + I2K  + 2I2}  2 (45) 

and  X = {(1  + ai)/I2K){(l-f)X2  + (1-f  +2I2K)X  + (l-f)X  (aj/b^/6  + 

2 

+ I2K+2I2}  / [[6  (1-f)  xWk  + 12X  + 12f/K  + 

+ (ai/bi)Xo(2(l-f)X/l2K  + 2}]  {(1-f)  X2  + (1-f  + 2I2K)X  + 

+ (1-f ) Xq  (ai/bi)/6  + + 2l2}  - [2(l-f)X3/^K  + 6X2  + 

+ 12f  X/K  + 12f/K2  + X (a  ,/b.)  {(l-fJX^^K  + 2X  + 2/K}  ] {2  (1-f ) X + 

0 l l 


+ 1 - f + 2I2K}]. 


Thus,  Q (0, T)  and  M(0,T)  may  be  evaluated  from  equations  (22) 
and  (23)  with  the  aid  of  equations  (41)  , (43)  , (45)  and  (46) . 

Hence,  Q(a,T)  and  M(a,T)  may  be  calculated  from  equations  (15), 
(16) , (19)  and  (20)  and  the  migration  of  the  generalized 
stresses  followed  along  the  facet  i during  this  phase  of  motion, 
which  is  completed  at  a dimensionless  time  Tk+1  , when  the 
generalized  stresses  reach  the  next  node  i+1  or  return  to  node  i. 
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5.  Impulsive  Loading  of  a Simply  Supported  Beam 

5.1  Introduction 

Consider  the  simply  supported  rigid  perfectly  plastic  beam 
shown  in  Figure  6(a)  which  is  subjected  to  a uniformly  distri- 
buted impulsive  velocity  V . It  was  found  in  Reference  [1] 
that  the  behavior  of  this  beam  is  governed  by  transverse  shear 
effects  alone  and  is  independent  of  rotatory  inertia  when 
0 <.  v <.  1 and  where 

v = QoL/2Mq.  (47) 

Furthermore,  it  may  be  shown  that  the  associated  generalised 
stresses  M and  Q throughout  the  entire  beam  either  lie  in  or 
within  the  Ilyushin-Shapiro  yield  curve.  Thus,  it  is  only 
required  to  examine  the  case  v M in  this  section.  A numer- 
ical solution  is  sought  using  the  Ilyushin-Shapiro  yield  curve, 
which  is  given  by  equations  (3)  and  illustrated  in  Fi' -re  2. 

The  response  of  this  beam  consists  of  two  phases,  which  are 
described  in  the  two  following  sections. 

5.2  First  Phase  of  Motion 

It  is  assumed  that  the  beam  in  Figure  6(a)  responds  with 
the  velocity  field  illustrated  in  Figure  6(b).  The  regions 
0 £ x — Zjl  and  z^  — x — L are  rigid,  while  a plastic  zone  de- 
velops within  the  region  z^  — x — z^. 

The  transverse  shear  force  (Q)  in  the  plastic  zone 
(z  £ x £ L)  is  quite  small  and  positive  in  the  theoretical 
analysis  presented  in  Reference  [1] . This  suggests  that  a 
linearisation  similar  to  equation  (6)  may  be  used,  or 
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M/Mq  = -1  + eQ/Qo  (48) 

when  z^  x i z^.  Thus,  equations  (10)  to  (12)  are  again  ob- 
tained which  give  equation  (13a)  and 

w(x,t)  = n ^ ( E - fF) coshDx  + ^ (F-fE) sinhftx  + VQ  (49) 

for  ^ x — when  satisfying  <Jj(x,0)  = 0 and  w(x,0)  = VQ 
and  where  ft  and  f are  defined  by  equations  (9b)  and  (17f) , 
respectively . 

The  velocity  fields  in  the  rigid  regions  0 £ x — z^  and 
z2  — x — L are 

w(x,t)  = C ( t)  + D(t)x  (50a) 

and  w(x,t)  = G(t),  (50b) 

respectively. 

It  may  be  shown  that  equations  (13a),  (49),  and  (50)  give 

• 

W(a,T)  = V,  + (V  - Vx) a (51) 

in  the  rigid  zone  0 — a — 1, 

W(a,T)  = 1 + (V  -1)  [cosh{K(  82-a61)  } + f sinh  (k  ( ) > 1 /B  , ( 52a) 

and  ij>(a,T)  = K ( 1-V)  sinh  (k  ( B2-<*6^)  }/B  , (52b) 

where  B = cosh  (k  ( S2~  B^)  ) + f sinh  (K  ( B2~  B^)  ) (52c) 

in  the  plastic  zone  1 — ot  1 B an<^ 

W(a,T)  = 1 + (V-l)/B  (53) 

for  the  central  rigid  zone  $2/ ~ a ~ when  satisfying 

w(0,t)  = V1#  w(zlft)  = V,  [w(z^,t)]  = 0,  [w(z2,t)l  = 0,  and 

( vp ( z 2 , t)  1 = 0,  and  where 
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a = x/zlf  $1  = Zl/L,  &2  = z2/L,  K = «L,  I2  = Ir/mL2, 

• • 

V = V/Vo'  V1  = Vl/Vo'  W(a,t)  = w ( x , t ) / Vq , jj7(a,T)  = Li|;  (x , t) /V  , 

and  T = Mot/mVoL2.  (54a-j) 

• • 

Now,  if  = Zj  = 0 during  the  first  phase  of  motion  (i.e., 
no  travelling  interfaces),  then  equations  (1)  and  (51)  for 
0 < a < 1 give 

— 2 — — — 

Q(a,T)  = a B^V  - + ctgjV./^v  + 1 (55a) 

and  M(a,T)  = a^2^  - V)/6  - - 2ag1v  + aI2(V-V1)  (55b) 

since  M(0,T)  = 0 for  simple  supports  and  Q(0,T)  = 1 for  plastic 
flow  to  develop  at  the  supports,  and  where 

• • 

Q(a,T)  = Q(x,t)/Qo,  M ( a , t)  = M(x,t)/MQ,  V(a,t)  = mL2V (x , t) /Mq , 

• • 

and  V^ct,!1)  = mL2Vx  ( x , t)  /Mq . (56a-d) 


Similarly,  for  the  plastic  zone  (1  1 a 1 g^g^)  • 

Q(a,T)  = -V(sinh{K(62-a61)}  + fcosh{K(82  -a81)}]/2vBK  (57a) 

and  M(a,T)  = -1  - fV  [sinh(K  ( } + fcosh(K(B2  -ag^  } ] /BK2  , 

(57b) 

where  the  constants  of  integration  were  found  by  demanding 
that  equation  (48)  was  satisfied  throughout  the  plastic  zone. 
Finally, 


Q(a,T)  = 


8 


2 


f/K)/2vB 


(58a) 


f ■ ■ — ■ 1 
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and  M(a ,T)  = -1  - V[a2&12/2  - - afg^/K  + 1/2  + 

+(f/K)2/2]/B  (58b) 

in  the  central  rigid  zone  E^^l  ~ a ~ 1/BX  when  [Q ( B2/81 ,T) ] = 0 
and  [M(62/g1 , T) ] = 0. 

Thus,  equation  (58a)  satisfies  the  symmetry  requirement 
Q(1/8X»  T)  = 0 when 

B2  = 1 - f/K.  (59) 

Moreover,  9w(x,t)/3x  at  x = (i.e.,  lMz^t))  in  the  rigid 

region  equals  iJj  ( z ^ , t ) in  the  plastic  zone  when 

• • • 

V - Vx  = -gjKV  s i nh { K ( B 2 " Bx)  }/B  . (60) 

Furthermore,  the  conditions  [Q ( 1 , T ) ] = 0 and  [M(1,T)]  = 0 may 
be  solved  to  predict 

V = (vBx2/3  - Bx/2  - 2I2v)  (Bx3/12  + BXI2  + 8xfD/2K2  + 

B12D/3K  + I2D/K)-1  (61a) 

and  Vx  = - ( 2v8^2/3  - 8^2  + 2vBxD/K  - D/K  + 2I2v  + 

+2vfD/K2) (B13/12  + I2B1  + fB1D/2K2  + DBX2/3K  + I2D/K)-1,  (61b) 

where 

D = (sinh{K(B2  -B^}  + fcosh{K(B2  - 6X)  } ] /B , (61c) 

and  8X  and  82  are  related  to  the  length  of  the  plastic  zone 


and  satisfy  equation  (60)  which  can  be  written  with  the  aid 
of  equations  (61)  in  the  form 
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k(v012/3  - Bj^/2  - 2I2v)  sinh  (k(B2  - 8-^  }/B  = - vB-^  - 

- 2v81D/K  + D/K  - 2vfO/K2)/B1.  (62) 

Now,  equations  (61)  predict 


V = 1 + VT  and  = 1 + V-^T 


(63a, b) 


since  V - - 1 at  T = 0.  The  first  phase  of  motion  is 

completed  when  = 0,  or 

Tj_  = -1/Vlf  (64) 

where  is  given  by  equation  (61b) , and  the  corresponding 
transverse  displacement  profile  and  angular  deformations  can 
be  evaluated  analytically  in  a straightforward  manner. 

5.3  Second  Phase  of  Motion 

The  second  phase  of  motion  with  ^ 0 and  z2  = 0 
commences  at  T = T^,  when  plastic  shearing  ceases  at  the 
supports  ( = 0).  Thus,  equation  (51)  for  the  rigid  zone 
0 — a — 1 becomes 


W(a,  T)  = Va,  (65) 

while  equations  (52a)  and  (52b)  remain  valid  for  the  plastic 
zone  with  1 — a — 82/B^,  an(*  ec3uat;*-on  (53)  remains  valid 

for  the  central  rigid  region  - a - V3,.  The  require- 

ment  that  9w(x,t)/3x  at  x = z^  in  the  rigid  region  equals 
>Mx,t)  in  the  plastic  zone  at  x = z^  gives 
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V = BXK  tanh{K(B2  - ) > [ 1 + (f+BjK)  tanh{K(B2  - 6X)  > J . 

Equations  (1),  (52a),  (52b),  (53),  and  (65)  together 
with  equation  (48),  the  simply  supported  boundary  cor  ions, 
lQ( a,T) ] = [M( a»T) ] = 0 at  a = B2/B1  and  [Q(cs,T)]  = 0 * , a =1 


predict 
Q ( os , T ) 

M(a,T)  = 


• • • • 

= (VBj^  - VBX)  (1  - a2)  /4v  - VD/2vK  - (V  - l)D2$1/2v 

(BjV  - B1V)(a3B1/6  - 8^/2  - aI2/B1)  + 

— _ 2 

+ VaDBx/K  + (V  - 1)  D (67b) 


for  0 — a — 1, 

• • 

Q(a,T)  = -{V/K  + (V  - l)DB1Hsinh{K(B2  - afij.)  } + 

+ f cosh  {K(B2  - aB1))l/2vB  (68a) 

• • 

M(a,T)  = -1  - {V  + (V-l)  Bj^DK}  (cosh{K(82  - aB-^  }/BK2  + 

+ fsinh(K(82  - aBx) }/BK2  - I2  cosh{K(82  - aB^ }/B ] (68b) 

for  1 <.  a <.  B2/B^/  and 

• • 

Q(a,T)  = (V  + (V-l)  KDB-^ ) ( osB^  - B2  - f/K) /2vB  (69a) 

• • 

M(a,T)  = (V  + (V-DKDB^tl2  - 1/K2  - (a^2  - 2aB1B2  + B22)/2 
-f  (B2  - csB^/Ki/B  - 1 (69b) 

when  82/B^  — a — 1/8^,  where 


mLV  z./M  . 
01  o 


(66) 


(67a) 


(70) 
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The  symmetry  requirement  Q(1/3^,T)  = 0 again  gives 
equation  (59),  while  [M(1,T)]  = 0 with  equation  (66)  predicts 

V[312/3  + DSj/K  + I2  + 1/K2  - I2  cosh{K(32  - 0 ) }/B]  - 

- (0.^3  + I2/31  + D2B  cosech{K(32  - 0^  }/K  + 

+ DB  cosech  {K(B2  - Bj^)  }/K261  - I2D  coth{K(02  - = -1 


which  again  using  equation  (66)  becomes 


(71) 


&1  = -Q2[[312/3  + DB1/K  + I2  + 1/K2  - I2  cosh{K(B2  - B1)}/B]P  - 

- [B12K  tanh{K(B2  - B1) }/3  + I2K  tanh{K(B2  - 0^  } + 

+ D2B1B  sech{K(B2  - Bx) } + DB  sech{K(B2  - 0^) )/K  - I2DK]qJ-1, 


(72) 


where 


P = K tanh{K(B2  - 0^ } + fKtanh  {K(B2  - 0^  } - 
-SjK2  sech2  (K(B2  - 6X) } 


and 


Q = 1 + (f  + BjK)  tanh{K(32  - 0^  } . 


(73a) 


(73b) 


Motion  finally  ceases  at  T = T2  when  V = 0,  or  0^  = B2 
according  to  equation  (66),  where  B2  is  given  by  equation  (59). 

Equation  (72)  must  be  integrated  numerically  to  give  0^ 
as  a function  of  T during  the  second  phase  of  motion.  The 
velocity  V and  acceleration  V may  then  be  calculated  from 
equations  (66)  and  (71),  respectively.  Hence,  the  dimensionless 
bending  moments  and  transverse  shear  forces  may  be  determined 
from  equations  (67)  to  (69) . The  total  transverse  displacement 
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profile  at  any  position  and  time  is  found  by  a straightforward 
integration  of  the  velocity  field  given  by  equations  (51), 

(52a),  (53),  and  (65)  when  making  due  allowance  for  a travelling 
interface  during  the  second  phase  of  motion.  In  particular, 
the  total  transverse  displacement  at  the  supports  is  accumu- 
lated during  the  first  phase  of  motion  since  no  further  plastic 
flow  develops  at  the  supports  during  the  second  phase  of 
motion,  i . e. , 


rTi  - 

W ( 0 ,T)  =J  VjdT  = 1^/2  , 


(74) 


where  and  are  given  by  equations  (63b)  and  (64) , 


T > T^,  and 


2t  2 


W ( 0 ,T)  = MQw(o,t)/mVo  L . 


(75) 


The  transverse  displacement  accumulated  at  the  mid-span  during 
the  second  phase  is 

• T, 


W(l/ei,T2)  - Wll/^,^)  = J 2 Wd/e^TjdT  , (76) 


where  W(1/$^,T)  is  given  by  equation  (52a),  and 


Wd/B^T^  = Tj.(l  + VT1/2B) 


(77) 


is  the  contribution  from  equation  (53)  for  the  first  phase. 
The  total  slope  at  the  supports  at  the  end  of  the  first  phase 
of  motion  is 

»T , 


§(0,T1)  = J 1 


(V  - V1)dT/B1  = (V  - V1)H1  /2Q1  (78) 


according  to  equation  (51),  while  the  additional  slope  accumu- 
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lated  during  the  second  phase  of  motion  is 

T 

6(0, T2)  = J 2 VdT/6x  (79) 

T1 

according  to  equation  (65),  where 

0(0, T)  = Mq6(  0,t)/mLVQ2.  (80) 

It  may  be  shown  that  the  various  equations  in  section  4 
of  Reference  [1]  for  the  problem  illustrated  in  Figure  6a 
with  a square  yield  curve  are  recovered  from  the  corresponding 
equations  with  e = 0 in  sections  5.2  and  5.3  here. 
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6.  Discussion 
6.1  Impact  of  a Mass  on  a Long  Beam 

The  numerical  predictions  for  the  response  of  the  rigid 
perfectly  plastic  beam  in  Figure  4(a)  hit  by  a mass  travelling 
with  an  initial  velocity  VQ  are  given  in  Figures  7 to  12. 

The  numerical  values  in  Figure  7 indicate  a relatively 
small  difference  between  the  results  for  2 facets  and  50  facets 
per  quadrant,  while  the  results  for  10  facets,  which  are  not 
shown  in  Figure  7,  lie  very  close  to  the  results  for  50  facets. 
The  theoretical  results  in  Figures  7 and  8,  with  transverse  shear 
retained  according  to  the  Ilyushin-Shapiro  yield  condition  but 
rotatory  inertia  neglected  (I  = 0),  are  compared  with  the  theo- 
retical predictions  of  Symonds  [3]  for  a square  yield  condition 
and  the  simple  bending  only  solution.  It  is  clear  that  trans- 
verse shear  effects  lead  to  a dramatic  reduction  in  the  angle 
underneath  the  striker  (|"0|)  and  a significant  increase  in  the 
maximum  transverse  displacement  underneath  the  striker  (W) . 
Moreover,  the  angle  ("0)  is  quite  sensitive  to  the  actual  shape 
of  the  yield  curve,  while  the  maximum  displacement  (W)  is  less 
sensitive. 

The  numerical  results  for  the  combined  influence  of  trans- 
verse shear  and  rotatory  inertia  on  the  dynamic  response  of  beams 
having  various  cross-sections  are  compared  in  Figure  9 with  the 
corresponding  theoretical  results  presented  in  Reference  [1]  for 
a square  yield  curve.  The  curves  labelled  (l)  to  (5)  in  Figure 


9 respectively  correspond  to  Wl4x87,  Wl2x40,  W24x55  wide  flanged 
I-sections  and  a 3/8  in.  wide  x 1 in.  deep  rectangular  cross- 
section  beam  hit  by  masses  with  mH/G  = 1.10651.  This  mass  ratio 
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was  used  by  Parkes  [11]  in  some  of  his  cantilever  tests.  H 
was  taken  as  the  total  depth  of  a beam,  and  Qq  was  calculated 
using  an  area  equal  to  the  beam  depth  times  the  web  thickness 
and  a shear  yield  stress  of  oq/2  as  discussed  in  Reference  [1] . 
Again  the  shape  of  the  yield  curve  exercises  an  important 
effect  on  the  angle  0,  while  the  maximum  transverse  displace- 
ment (W)  is  fairly  insensitive,  particularly  for  the  larger 
values  of  \Q.  Thus,  Figures  8 and  9 suggest  that  an  analysis 
with  a square  yield  curve  (Reference  [1]  when  I/O  and  Reference 
[3]  when  1=0)  is  adequate  when  the  maximum  transverse  displace- 
ments are  of  interest.  However,  the  Ilyushin-Shapiro  yield 
curve  must  be  used  when  the  rotation  under  the  striker  (0)  is 
required. 

The  theoretical  results  in  Figure  10  indicate  the  relative 

importance  of  retaining  transverse  shear  forces  in  the  yield 

criterion  and  rotatory  inertia  in  the  basic  equations.  The 

influence  of  transverse  shear  forces  is  consistent  with  the 

earlier  remarks,  while  rotatory  inertia  is  not  noticeable  in 

the  curves  for  XQ  = 24.357  (e.g.,  Wl4x87  wide  flanged  I-section 

with  mH/G  = 1.10651),  but  does  lead  to  a reduction  in  0”  of 

10%,  approximately,  when  X = 3.31953  (e.g.,  3/8  in.  wide  x 

o 

1 in.  deep  rectangular  beam  with  mH/G  = 1.10651). 

The  partition  of  the  total  energy  and  migration  of  the 
generalised  stresses  around  the  Ilyushin-Shapiro  yield  curve 
are  shown  in  Figures  11  and  12,  respectively,  for  XQ  = 8.102 


and  I = 0.4276. 
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6.2  Impulsive  Loading  of  a Simply  Supported  Beam 

The  numerical  predictions  for  the  dynamic  plastic  response 
of  a simply  supported  beam  governed  by  the  Ilyushin-Shapiro 
yield  curve  with  e = 0.06  and  subjected  to  a uniform  impulsive 
velocity  Vq  as  shown  in  Figure  6(a)  are  presented  in  Figures 
13  to  17.  The  curves  labelled  (l)  in  Figures  14  and  16  corre- 
spond to  a beam  with  a rectangular  cross-section,  and  the 
curves  labelled  (2)  , (5)  , and  (4)  are  associated  with  the  wide- 
flanged  I-sections  having  the  cross-sections  respectively 
labelled  (3),  (2),  and  (T)  in  Figure  9.  The  results  in  Figures 
13  to  16  are  virtually  indistinguishable  from  the  corresponding 
theoretical  values  presented  in  Reference  [1]  for 
a square  yield  curve.  Thus,  the  simpler  theoreti- 
cal predictions  in  Reference  [1]  for  a square  yield  curve  are 
adequate  for  the  problem  shown  in  Figure  6(a).  The  presence 
of  a central  rigid  zone  — a — 1/3^  which  was  considered 

in  sections  5.2  and  5.3  does  not  appear  to  influence  the  results 
significantly.  Indeed  &2  = 0.985  for  the  particular  beam  ex- 
amined in  Figure  17  which  means  that  the  central  rigid  zone  has 
a total  width  of  only  0.03L  in  this  case.  The  duration  of  the 
first  phase  of  motion  is  T^  = 0.10765  which  is  to  be  compared 
with  T^  = 0.10795  found  in  Reference  [1]. 

It  may  be  shown  that  the  generalised  stress  profiles  associ- 
ated with  the  theoretical  solutions  having  I = 0 in  References 
[1]  and  (4]  for  a square  yield  curve,  lie  within  or  on  Robinson's 
circular  yield  curve  [5,10]  which  in  turn  inscribes  the  Ilyushin- 
Shapiro  yield  curve.  Thus,  the  theoretical  results  with  1=0 
in  section  4.1  of  Reference  [4]  and  in  Reference  [1]  are  identi- 
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cal  to  the  corresponding  theoretical  predictions  using  an 
Ilyushin-Shapiro  yield  curve,  while  the  theoretical  predictions 
in  Reference  [1]  with  I/O  are  virtually  indistinguishable 
from  the  numerical  results  obtained  herein. 

It  is  evident  from  Figure  14  for  v = 1.5  that  the  retention 
of  transverse  shear  effects  in  an  analysis  with  1=0  gives 
the  same  values  of  W (dimensionless  transverse  displacement  at 
mid-span)  which  are  predicted  by  a simple  bending  theory,  while 
the  consideration  of  I leads  to  a reduction  in  the  mid-span 
transverse  displacements  (W)  up  to  ten  per  cent,  approximately. 

On  the  other  hand,  the  inclusion  of  transverse  shear  effects 
causes  a significant  reduction  in  0 (dimensionless  angle  of  ro- 
tation at  supports),  while  consideration  of  rotatory  inertia  (I) 
is  responsible  for  a further  decrease  up  to  17  per  cent,  approx- 
imately. 

The  energy  partition  in  Figure  15  for  v = 1.5  is  changed 
significantly  when  the  influence  of  transverse  shear  effects 
is  included  in  the  yield  curve  since  fifty  per  cent  of  the  ini- 
tial energy  is  then  absorbed  due  to  shearing  deformations  alone 
when  T 2.  1/6.  The  incorporation  of  rotatory  inertia  (I)  in  the 
analysis  causes  a further  increase  in  the  amount  of  energy  ab- 
sorbed due  to  transverse  shear  forces. 

It  is  evident  from  the  numerical  results  in  Figure  16  that 
transverse  shear  effects  and  rotatory  inertia  have  a negligible 
influence  on  the  maximum  permanent  transverse  displacements  at 
the  mid-span  when  v>4,  approximately. 
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6.3  General  Comments 

In  order  to  avoid  complete  severance  of  the  beams  at  the 
shear  hinges  it  is  necessary  to  ensure  that  the  maximum  shear 
displacements  in  the  numerical  studies  are  less  than  some  pro- 
portion of  the  beam  thickness  [12] . 

The  material  was  idealised  as  rigid  perfectly  plastic  in 
this  article  since  this  characterisation  gave  satisfactory 
agreement  between  theoretical  predictions  and  experimental  re- 
sults for  many  problems  when  the  external  dynamic  energy  was 
greater  than  about  five  times  the  strain  energy  which  corld  be 
absorbed  in  a wholly  elastic  manner  [13] . In  addition,  the 
influences  of  material  strain  rate  sensitivity  and  geometry  changes/ 
or  finite  transverse  displacements,  were  disregarded.  These 
assumptions  are  probably  reasonable  for  impulsively  loaded 
strain  rate  insensitive  beams  when  supported  without  axial  re- 
straints, and  for  infinitely  long  strain  rate  insensitive  beams 
struck  by  a mass  when  the  maximum  transverse  displacements  are 
less  than  the  beam  thickness,  approximately  [13].  The  influence 
of  material  strain  hardening  has  been  neglected  in  this  work, 
but  some  general  remarks  on  this  topic  are  presented  in  Refer- 
ences [13 , 14]  . 

The  numerical  procedure  developed  in  this  article  could 
be  used  to  examine  the  behavior  of  the  two  beam  problems  when 
governed  by  any  convex  yield  criterion  which  relates  the  trans- 
verse shear  force  and  bending  moment  required  for  plastic  flow. 
Moreover,  the  numerical  scheme  could  be  employed  to  solve  many 
other  beam  problems  in  an  inexpensive  manner  and  appears  suffi- 
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ciently  attractive  to  warrant  further  development  in  order  to 
study  the  dynamic  plastic  response  of  plates  and  shells. 


7.  Conclusions 

A numerical  procedure  has  been  used  to  examine  the  in- 
fluence of  transverse  shear  forces  in  the  yield  criterion 
and  rotatory  inertia  on  the  dynamic  plastic  response  of  beams. 
Various  results  are  presented  for  a long  beam  impacted  by  a 
mass  and  a simply  supported  beam  loaded  impulsively,  both  of 
which  are  made  from  a rigid  perfectly  plastic  material  with 
yielding  controlled  by  the  Ilyushin-Shapiro  yield  criterion. 

Transverse  shear  effects  lead  to  a dramatic  reduction  in 
the  slopes  of  the  deformed  profiles  for  both  beam  problems. 
Moreover,  the  slope  of  the  deformed  profile  underneath  the 
striker  in  the  impact  problem  is  quite  sensitive  to  the  actual 
shape  of  a yield  curve,  while  the  maximum  transverse  displace- 
ment is  less  sensitive.  The  retention  of  rotatory  inertia 
in  the  basic  equations  may  lead  to  a further  10  per  cent 
reduction  in  the  slope  of  the  impact  problem,  and  reductions 
of  up  to  17  per  cent  and  10  per  cent  for  the  slopes  and  trans- 
verse displacements  of  the  impulsive  problem,  respectively. 


, 
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I = 0.1924,  x indicates  T = T^ • 

Figure  16.  Variation  of  Maximum  Permanent  Transverse  Displace- 
ment (WT)  of  Simply  Supported  Beam  in  Figure  6(a)  with  v 
and  e = 0.06. 

• simple  bending  theory 

1 = 0 and  e = 0 (square  yield)  [4]  . 
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0 Iv  = 0.2886,  @ Iv  = 0.1584, 

0 Iv  = 0.0781,  0 Iv  = 0.0600 

Figure  17.  Migration  of  Generalised  Stresses  in  the  Simply 
Supported  Beam  in  Figure  6(a). 
v = 2,  I = 0.1443,  e = 0.06,  = 0.10765. 
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A numerical  procedure  is  used  to  examine  the  influence  of 
transverse  shear  forces  in  the  yield  criterion  and  rotatory 
inertia  on  the  dynamic  plastic  response  of  beams.  Various 
results  are  presented  for  a long  beam  impacted  by  a mass  and 
a simply  supported  beam  loaded  impulsively,  both  of  which  are 
made  from  a rigid  perfectly  plastic  material  with  yielding 
controlled  by  the  Ilyushin-Shapiro  yield  criterion.  
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and  10  per  cent  in  the  slopes  and  maximum  transverse  displace- 
ments, respectively. 
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